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The equilibrium properties of flux lines pinned by columnar disorder are studied, using the 
analogy with the time evolution of a diffusing scalar density in a randomly amplifying medium. 
Near Hci , the physical features of the vortices in the localized phase are shown to be determined by 
the density of states near the band edge. As a result, Hci is inversely proportional to the logarithm 
of the sample size, and the screening length of the perpendicular magnetic field decreases with 
temperature. For large tilt the extended ground state turns out to wander in the plane perpendicular 
to the defects with exponents corresponding to a directed polymer in a random medium, and the 
energy difference between two competing metastable states in this case is extensive. The divergence 
of the effective potential associated with strong pinning centers as the tilt approaches its critical 
value is discussed as well. 



PACS numbers: 74.60. Ge, 72.15.Rn, 05.70.Ln 
I. INTRODUCTION 

The physical properties of vortex hnes in the mixed 
phase of type II superconductors has become subject of 
intense research in recent years When an external 
current density is apphed to the bulk of the system, the 
flux lines may start to move under the action of the 
Lorentz force. Within a perfectly homogeneous system 
this driving Lorentz force is balanced only by the friction 
force opposing the steady state velocity of the flux lines, 
so that there is a dissipation coupled to the appearance of 
a finite electric field which results from the fiux motion. 
To avoid this effect, flux lines in the mixed phase should 
be pinned by inhomogeneities in the underlying crystal 
structure, such as the point defects associated with va- 
cancies of oxygen atoms in high temperature cuprates 
. It turns out that the pinning of the vortices is much 
stronger when these impurities are in the form of corre- 
lated disorder, like twin boundaries or columnar defects 
, aligned along the direction of the external magnetic 
field. If the thermal fluctuations are small enough, the 
flux line may lie along the entire extended defect in the 
bulk; this is in contrast to short scale disorder pinning, in 
which the line should accommodate itself to the potential 
fluctuations and hence increases its elastic energy. How- 
ever, the correlated defects pinning becomes less effective 
when the direction of the external magnetic field is tilted 
with respect to the anisotropy axis, which we take to be 
along the z direction. At some critical tilt, for which 
the free energy per unit length of the defect is less than 
the free energy associated with the perpendicular field, 
a pinning-depinning phase transition occurs and the flux 
lines are delocalized. 

The static and dynamic response of the flux lines in 
the presence of columnar defects have been considered 
by Nelson and Vinokur Using the mapping of flux 
lines in d -t- 1 dimensional superconductor to the world 
lines of bosons in a d— dimensional quantum system, the 
authors identifled the phase space diagram of the system 
which contains a high temperature "superfluid" and low- 



temperature "bose glass" phases, as well as Mott insula- 
tor at the matching field, = rtpm^Oj for which there is 
one fiux line per defect. At low temperature, this match- 
ing field separates the "dilute" region of the bose glass 
phase, for which the vortex lines are pinned individually 
by the defects, from the high density region, where in- 
teractions are important in determining the localization 
length and transport properties of the flux lines. 

In the low field region, the vortices are localized by the 
interaction with the correlated defects. Each pin is the 
analog of a 2D potential well which may be described (up 
to logarithmic corrections) as a cylindrical square well. 
The temperature of the superconductor, in turn, corre- 
sponds to the Plank constant h of the quantum boson 
system. For the dilute vortex arrays, where the pinning 
energy is larger than the interaction energy, there are two 
regimes. For low temperature the localization length is 
given approximately by the radius of the defect so that 
each flux line is localized by one defect. As the tempera- 
ture increases, the localization length of one defect grows 
exponentially with T^, and the flux line is then localized 
by several defects, forming an effective d dimensional po- 
tential well in the corresponding boson system. 

The depinning of flux line which occurs as a result 
of external field tilt has been carefully investigated by 
Hatano and Nelson The Hamiltonian of the cor- 

responding boson problem is no longer Hermitian; the 
kinetic term of the Hamiltonian is subject to an imag- 
inary gauge transformation where the gauge field h is 
related to the perpendicular magnetic field Hj^ via h = 
H^(/)o/(47r), where ipQ is the fiux quantum. In the ab- 
sence of tilt, the probability of finding a point of the 
vortex at transverse displacement r relative to the center 
of the pin is independent of z, and given by the square 
of 0gs(r), the ground-state wavefunction of the Hermi- 
tian Hamiltonian. For small tilt, there are left and right 
ground-state eigenfunctions which correspond to the left 
and right "tilting" of the Hermitian ground-state are still 
localized. It turns out that the flux line is described by 
these left and right eigenfunction at the bottom and the 
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top of the sample, while deep in the bulk the probabil- 
ity function approaches its Hermitian limit. Typically, 
the "surface roughness" associated with the tilt extends 
into the bulk up to some characteristic distance which 
diverges as the tilting angle approaches the critical an- 
gle, for which the flux line delocalizes and the current 
response becomes linear. 

This non-Hermitian delocalization process has been 
shown to be of importance for other physical systems 
as well. In general, the time evolution of diffusing scalar 
field in random environment is determined by the eigen- 
states of Liouville operator which is the analog of the 
Hermitian Hamiltonian of a quantum particle in a disor- 
dered system; the effect of convection may be modeled by 
the imaginary gauge transformation. The appearance of 
extended states as convection is increased beyond a crit- 
ical value, associated with the complex spectral points 
of the resulting non-Hermitian operator, is the manifes- 
tation of the delocalization transition. The dynamics of 
strongly driven charge density waves , and the growth 
modes of biological populations are among the sys- 
tems considered. 

In this paper we study this tilt induced delocalization 
transition at very dilute flux line concentrations, i.e., near 
Hci- We assume, therefore, that only the low- lying free 
energy states of the system are occupied, and that the 
repulsive interaction between vortices may be taken into 
account by sequential filling of energy levels. Although 
this assumption is self consistent in the localized phase, 
it turns out that it is valid in the extanded phase only if 
the sample is finite. If the size of the sample is taken to 
infinity at finite vortex density, interaction may lead to a 
smecticlike flux crystal, as has been discussed in Ref. [|j . 

This paper is organized as follows. In Section 2, a brief 
introduction to the boson - vortex analogy is presented, 
as well as the correspondence between the statistical me- 
chanics of flux lines and the time evolution of a scalar 
field. Section 3 deals with the finite size effects in the 
localized phase; both the critical field Hd and the pene- 
tration depth of the perpendicular field are shown to be 
determined by the statistics of the tail of the density of 
states. Since this function is strongly related to the size 
of the system, physical properties of the pinned vortex 
are determined by the actual size of the bulk. In Section 
4 we consider the delocalized phase, in which we recog- 
nize that the low lying states are extended and are related 
to the properties of Burgers's equation with conservative 
noise in 1 + 1 dimension ||^. Some technical details are 
given in the appendices. 



II. FLUX LINE IN COLUMNAR DISORDER 



We review here the basic physics of a single flux line in 
d -|- 1 dimensional superconductor, as has been discussed 
in |] and §. 




FIG. 1. Vortex line in a superconductor with columnar 
disorder. If H±_ = 0, the flux line is localized, i.e., trapped 
by one or few pins into some region in the plane perpendic- 
ular to the correlated disorder. As the external magnetic 
field is tilted away from the columns, the flux line tends 
to delocalize and tilt in the direction of the external field. 



Consider a flux line in a superconducting sample of 
thickness Lq when the sample is pierced by columnar 
pins, such as long alighted columns of damaged material, 
as illustrated in Fig. 1. 

Neglecting overhangs of the line, we define the flux line 
by its trajectory r(z). The free energy of this line may 
be written as. 
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where ei is the tilt modulus of the flux line, i.e., the 
energy of the flux line per unit length, and the clas- 
tic contribution ^(^^^^)^ is the flrst nontrivial term in 
the small tipping angle expansion of the line energy of a 
nearly straight vortex line. U{r) is the random potential 
which arises from a ^-independent set of disorder-induced 
columnar pinning potentials, and Hj^ is the perpendicu- 
lar field. U{r) is taken from some bounded distribution, 
such as square distribution of width A around zero; its 
bias, a, is determined by ei and the parallel component 
of the external magnetic field, i.e., a = (ei — -j^)- Note 
that we assume H± to be much smaller than H|| , hence 
we neglect terms of order in (2.1). 

Using the analogy between the partition function 
for this flux line, Z = e~^/-^, and the path integral 
formulation of quantum mechanics, it is easy to see 
that the conditional probability Z{r, z;ro, zq) obeys the 
Schroedinger-like equation. 
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U{v)]Z -^h.^-VZ (2.2) 
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is the dimensionless perpendicular 



Let us now map this flux line model, which satisfies Eq. 



(2.2), onto a model for the diffusion of a passive scalar 
density, 4>{x), suppressed/ amplified by a quenched ran- 
dom potential, and driven by an external drift. The time 
evolution of such a system is described by the equation: 



dtcj)iv, t) = 7?V20(r, t)^[b+ l/(r)](/.(r, t) + h • V(; 



(2.3) 



where V{v) is the random amplification/suppression rate, 
and D is the difFusivity. Clearly, the partition func- 
tion of the flux line satisfies Eq. ( |2.3| ) for </>, where D 
plays the role of temperature divided by the tilt modu- 
lus, D — T/(2ei), h is proportional to the perpendicular 
magnetic field, h = ^ and the columnar disorder in 
the superconductor, as well as its bias, are normalized by 
the temperature, i.e., V'(r) = U{r)/T and b — a/t. The 
correspondence between the various scalar field and mag- 
netic vortex quantities are summarized in table 1. 

Table 1 
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Flux Line 
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Eq. ( |2.3| ) may be written as dt4> = —Ticf), where Ti., the 
Hamiltonian, is the linear operator which generates the 
time evolution of (f>, 



[b+V{r)] - ii - V. 



(2.4) 



The dynamics of this system is determined by the eigen- 
values and the eigenvectors of the Hamiltonian. For 
h = 0, the operator, H = -D'\/'^ + [b+V{r)], is Hermitian 
and it is well known that, for strong enough disorder, all 
its eigenf unctions are real and localized; the localization 
length is maximal at the center of the energy band, and 
minimal at the tails |10|. When h 7^ 0, the above Hamil- 
tonian is no longer Hermitian, and it may be diagonlized 
using a system of left and right eigenvectors. Since the 
perpendicular field term may be absorbed into the Lapla- 
cian by completing the square, V V -|- the right 
and left eigenfunctions of the new Hamiltonian are re- 
lated to the eigenfunctions of 7Y(h = 0) via an imaginary 
gauge transformation; if 0jj^''(r) is an eigenfunction of 
the Hermitian problem, then 
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are the eigenfunctions of the non-Hermitian operator 
with the same eigenvalue, provided that the local- 
ization length in the non-driven problem, is less than 
D/h. Thus, for small tilt, there is a spectral rigidity - 
the eigenvalue spectrum is robust. When h increases, 
the eigenfunctions become extended and the boundary 
conditions of the sample should be taken into account. 
Since (2.5), in general, does not satisfy these boundary 
conditions, the new eigenfunctions are no longer related 
to the h = case by a simple gauge transformation. 
As a result, the eigenvalues are changed; in case of pe- 
riodic boundary conditions, complex eigenvalues appear 
when D/h becomes smaller then As h is increased, 
these delocalized states appear first at the band center, 
for which the localization length is maximal, then move 
outwards, as has been discussed in P|. 



The mapping of Eq. (2.2) to Eq. (2.3) implies that for 
thick samples, Lq ~^ 00, the thermally excited flux line 
relaxes into the spatial configuration equivalent to the 
ground state eigenstate of Ti. For small h^, this state is 
localized, so that the flux line is localized around some 
spatial point at the bulk of the sample. This phenomenon 
is known as the transverse Meissner effect. As has been 
shown by the ends of the pinned flux line begin to 
tear away from the pinning center as the tilt increases. 
The "penetration depth" of the perpendicular field, as- 
sociated with the width of the region near the surface 
in which the transverse Meissner effect breaks down, di- 
verges as the tilt approaches he. Above this value, the 
eigenstate delocalizes, and the resulting trajectory wan- 
ders across the entire sample. 



III. MAGNETIC RESPONSE IN THE PINNED 
PHASE 



A. Size dependence of Hd 



n:h ' 



n,h=0 



The lower critical field, Hd, is defined in a random 
system as the minimal field in which the first vortex line 
enters the sample. One could easily recognize from Eqs. 
(2.1-2.4) that the condition for one flux line in the sample 
to be energetically favored over the completely diamag- 
netic phase, is that the Hamiltonian (2.4) admits nega- 
tive eigenvalues; the spectral points of (2.4) are propor- 
tional to the energies per unit length of the corresponding 
conflgurations of the vortex. Hd depends, thus, on the 
ground states of the Hamiltonian; for random potential 
one should consider the statistical properties of the lower 
tail of the density of states. Note that we consider here 
the physical case of bounded disorder, such that for each 
(2.5) sample there is a well defined lower critical field. This 
is in contrast to the discussion given by Larkin and Vi- 
nokur, who address the question of magnetic response 
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for unbounded, Gaussian random potential in an infinite 
sample. In such a case there is no well defined Hd [pl| . 

For a statistically homogeneous sample, any rare fluc- 
tuation is realized as one takes the bulk size to infinity. 
In particular, one may finds any large region in which 
the potential is close to its lower bound, so for infinite 
sample the minimal value of the potential V{x) sets the 
magnitude of Hd, i.e.. 
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For finite sample this is not the case. Typically, as the 
sample size increases, Hd — > Hd', however, the form of 
the size dependence of the critical magnetic field is not 
trivial. 

Consider, for example, a hole in the superconductor, 
i.e., a region of linear size L where the elastic energy of 
the flux line vanishes, ei = 0. The free energy of a vortex 
line which is trapped in this hole comes entirely from the 
limits on its transverse motion sets by the boundaries of 
this hole. Far away from the bulk critical temperature 
(i.e., when ii in the bulk is much bigger than the typical 
free energy per unit length associated with the trapped 
flux line) one may assume no vortex displacement along 
the hole boundaries, and the situation corresponds to a 
quantum mechanical particle in a box. In such a case the 
tail of the spectrum of (2.4) looks like 



e{k) - Dk^ 



(3.2) 



where the minimal value of k is sets by the system linear 
size L, kmin ~ ^/L, so that. 
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(3.3) 



This result reflects the fact that the entropy of the vortex 
line decreases if it is confined to finite region, hence its 
free energy per unit length increases. At zero tempera- 
ture this finite size correction to Hd disappears; as the 
temperature approaches its critical value, the potential 
steps at the hole boundaries become small, and in 2D this 
contribution to the free energy vanishes exponentially. 

Let us consider now the disordered case. In Appendix 
A we show that for a simple lattice model with square 
distribution of the disorder the density of states per unit 
volume near the minimal possible value of the energy 
(which we set, for convenience, to zero since this value is 
absorbed in the definition of Hd) is 



5(e) 
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(3.4) 



where Ig is the lattice constant. One may estimate the 
energy difference Ae between the ground state and zero 
for a typical sample of linear size L by the condition 



(3.5) 



so for large samples the ground state energy is given by 

D 



and 
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(3.7) 



In contrast to the hole example, one sees that the critical 
parallel field approaches Hd much slower, i.e., finite size 
effects are much stronger in the presence of disorder. For 
the tilt field below its critical value the situation remains 
the same, since the spectrum does not change. 

Upon taking reasonable parameters for high temper- 
ature superconductors, such as T ^ IQ^K, Iq ~ 100 — 

o o 

1000 A and A ~ 1000 A, one finds, however, that the 
entropic corrections to Hd are very small, 5Hd/Hd ~ 
10^^. Thus, for almost any realistic sample size one could 
take the value of Hd as the lower critical field. Since in 
most cases the radius of the columnar pins is comparable 
to, or bigger than, the superconducting coherence length, 
the range of the external field corresponding to the Meiss- 
ner phase turns out to be very narrow. For other physical 
situations when the time evolution of some scalar density 
is given by (2.3), these finite size corrections may by of 
importance. 



B. Transverse Meissner Effect for Dilute 
Concentration 

Given a system described by the equation 



(3.8) 



where Ti. is, in general, a time independent non-Hermitian 
operator with a complete set of left and right eigen- 
states (j)n,L,4'nM with the corresponding eigenvalues {e„} 
(i.e., the complex "energy spectrum" of the possibly non- 
Hermitian operator), the time evolution of the normal- 
ized amplitude with initial condition 0(r, 0) is given by 
its spectral decomposition: 



(r,t) 



(3.9) 



/ dr 0(r, t) 

^ En < '/'n,L |0(r,O) > (j)nnix)exp{-ent) 

S dvY,n< I '/'(r,0) > (?;'„„(a;)exp(-e„i) 

where < .. > is defined as the inner product < >= 
j dr tp*{r) (f>{x). At long times, the system is dominated 
by the ground state, i.e., the state for which the real part 
of e„. Re e„, is minimal. Using the mapping between the 
partition function of the vortex system and the scalar 
field (j), one observes that Eq. (^|^) describes the surface 
roughness of the fiux line: although deep in the bulk. 
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which is the equivalent to the "long time" of (3.9), the 
line is stuck to the spatial configuration corresponding 
to the ground state, this is not the case near the surface, 
where the "short time" limit of (3^) is relevant [||. 

Let us discuss this phenomenon for the h = case, 
when all the wavefunctions are localized. Spectral de- 
composition of small, localized, initial conditions near 
the most rapidly growing state requires projection upon 
the eigenvectors of the Hamiltonian. These projections, 
in general, decrease exponentially with the distance be- 
tween the fluctuation center and the localized eigenvec- 
tor. For example, if the initial fluctuation is centered 
around zero, the projection on an eigenstate localized at 
X is 



< (/<x=o(r,O)|0x(r) >~ exp-«"l^l, 



(3.10) 



where ac„ is the inverse localization length of the n-th 
eigenstate. The fastest growing states of TL dominate the 
system at long times; at any finite time, however, there 
is a competition between the initial conditions and the 
growth rates, which we explore using the following model. 

Assume for simplicity that the localization length is 
constant for all states of the system, i.e., it is indepen- 
dent of the eigenenergy as well as the spatial position. 
The probability amplitude (/)x generated by a localized 
initial condition '/'x=o(r,0) is proportional to 



0x(O 
</'x=o(i) 



(3.11) 



where ex is the eigenvalue of a growth mode localized at 
X. We shall study the relaxation of the first moment of 

(^(X,t), 

f - , N f dx X q!)(x, t) 

.(t).y.xx^(x,0.i^^^-^ (3.12) 

to the value x„s it assumes at long times, 



= / dx x (x) 



(3.13) 



where 0gs(x) is the lowest energy eige nmod e of the 
Hamiltonian. Given ex, one can optimize (3.11) in order 
to find the location of the eigenstate which dominates 
the system at finite time. For a finite sample it is clear 
that as i ^ oo the ground state (lower free energy) mode 
localized at position Xgg will dominate. The long time 
convergence of the normalized amp litude to the ground 
state location is also given by ( 3.11 ), where we use a nor- 
malization with respect to the final location amplitude, 
i.e., 



0x(t) 

'x = Xg4t) 



(3.14) 



In random samples one may assume that the eigenener- 
gies are described by a statistically homogeneous density 



of states 5(e), and a typical spatial distance R = |x— Xgs| 
associated with energy difference Ae = ex ^ — ex is given 
by the analog of Eq. (3.6) Q 



Ae R'' 5(e) = 1. 



(3.15) 



A similar mapping into variable range hopping in semi- 
conductors has already been exploited in Q for finite 
density of vortices in the superconducting medium. The 
long range repulsive interaction between flux lines has 
been taken into account, in the strongly localized regime, 
by forbidding multiple occupancy of localized states, and 
the "ground state" is defined by filling the energy band 
up to the "Fermi surface" determined by the chemical po- 
tential of the particles. If all the states are localized and 
the filling factor is such that the Fermi surface is not at 
the tails of the band, the density of states at the chemical 
potential, 5(/i), may be taken as a constant independent 
of the en ergy. Upon substituting Ae(i?) = l/{g{ii)R'^) 
into Eq. ( |3.14 ) , one can study the long time convergence 
to the ground state location, in terms of the distance 
R* away from Xg^ for which the normalized amplitude is 
maximal. One then finds that the approach of |x — Xg^l 
to zero at long times is given by a stretched exponential 



exp[ 



-(^) — ] 



(3.16) 



where r* is the characteristic time of the decay, t* — 
g{lJL)/K'^. Since the time in Eq. (2.3) corresponds to the 
the length in the correlated direction of the flux problem 
normalized by the temperature, the penetration depth of 
the perpendicular magnetic field is Tg{fi)/ k'^. 

For a very dilute concentration of flux lines, the ground 
state (to which the system eventually evolves), is in the 
lower tail of the density of states and is always localized, 
even for weak disorder in high dimensions. The nature 
of the surface roughness is determined by the other low 
energy states at the tail of the density of states func- 
tion. Since the ground state is now near the bottom of 
the band, the assumption that the density of states in 
this region is energy independent is no longer valid. One 
may estimate the time needed for the decay of a state 
in the tail to the gro und s tate by optimizing (3.14) with 
respect to R, using ( 3.15|) a nd (3.4). Taking the loga- 
rithm of both sides of ( 3^151 ), we get, up to logarithmic 
corrections, 
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(3.17) 



From this we obtain the semiclassical expression for 
the distance away from the ground state at time t: 



exp[ 



t/t* 



log (t/t*) 



1+2/d' 



(3.18) 



where now t* = l^/D, or in the vortex language, the 
width of the surface roughness is proportional to eiZg/T. 
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Note that the effect of temperature (or the diffusion con- 
stant) in this expression is to decrease the surface rough- 
ness rather then increasing it as in the former case; the 
reason for this counter-intuitive result is that the tem- 
perature dictates the energy spacing at the tail. Lower 
temperature implies bigger spacing, so that the effect of 
neighboring sites on a pinned flux line becomes smaller. 
Of course, the width of the tail its elf is proportional to 
the temperature, so that the result (3.18) becomes mean- 
ingless at the limit T ^ 0. 

This result implies that at small concentrations, i.e., 
close to Hci, the localization of flux lines is more effec- 
tive than in the case considered by Samples which 
are not thick enough to localize the vortices at high mag- 
netic field may do so near the lower critical field, when 
the effects of tail statistics are of importance. 

For a finite but small tilt h, the low-lying states are 
still localized, and the situation is essentially the same. 
At short times, there can be transition effects due to ex- 
tended states which overlap with the initial condition, 
but at long times the normalized amplitude is localized 
into the ground state. For strongly localized states near 
the band edge the tilt can be absorbed into the gauge 
transformation ( |2.5| ), the only change introduced in the 
consideration above is that the inverse localization length 
K becomes anisotropic, 



(h) = K(h 0) - h. 



(3.19) 



As a result, one should optimize ( 3.11 , 3.14 ) with respect 
to 9, the angle between h and x. The result turns out 



to be equivalent to replacing h in (|3.19| ) by its absolute 
value h. 



IV. EXTENDED STATES AT LARGE TILT 
A. Delocalized flux line as directed polymer 

Let us consider now the limit of very large tilt, in the 
sense that even the ground state of our system is delocal- 
ized, i.e., h » D/^Q. As ft- — > oo, the external magnetic 
field is perpendicular to the columnar defects; the flux 
line is now free in the z direction, while in the xy plane 
it sees the cross section of these defects, i.e., point dis- 
order. One might then guess that the flux line freely 
diffuses in the z direction, while in the xy plane it looks 
like a directed polymer in a 1 -I- 1 dimensional random 
medium, that is, the projection of the flux line on the xy 
plane is a pinned string characterized by the wandering 
exponent l/z = 2/3[^. It is separated from competing 
metastable states by energy barriers which (up to loga- 
rithmic corrections) scale like L^, where f3 — 1/3 and L 
is the linear system size [p|jl^] . 

The same features characterize the mo st rapidly grow- 
ing eigenmode of the Hamiltonian (2.4). The low ly- 
ing eigenstates of the non-Hermitian Hamiltonian are 
extended in the tilt direction while wandering in the 



transverse direction with the same exponents 1/z and 
(3. Suppose L is the length of the sample in the h di- 
rection, and W is the width in the perpendicular direc- 
tion, measured in units of some short length cutoff. For 
W » L^^^ there are many, spatially uncorrelated, com- 
peting metastable states in the sample; the decay of such 
a state into the ground state involves activation above 
a potential barrier of order L^/^ |]l3|, so that the time 
needed for such a process is of order io ''^ exp(L^/'^). On 
the other hand, as Vl^ < L^/'^, there is only one effective 
low energy state; since the energy barriers of the directed 
polymer are equivalent to the surface roughness in KPZ 
model, the potential barriers scale as W^^"^ in this case. 
In general, one sees that, in 2 -f 1 dimensions, when the 
low-lying states are extended, the time needed for the 
relaxation to the ground state is macroscopic and pro- 
portional to the (longitudinal or transverse) size of the 
system. 

In order to quantify these considerations, we may map 
Eq. (2.3) into a nonlinear problem with additive disorder 
using Cole-Hopf transformation, i.e.. 



(x, t) — exp^" 



(4.1) 



Separating the result into longitudinal terms, i.e., paral- 
lel to the drift h, and the transverse terms, perpendicular 
to h, one gets. 



dtp{x,Y^,t) = D,dlp{x.Y^,t) ^ D^Vlp{x,Y^,t) (4.2) 

+ - [Vip(2;, r^,t)Y + - [d,p{x, rx , t)Y 
7 7 

- \i - dxp{x,v^,t) + ^{b + V{x,Vperp)). 



Where for concreteness we take h along the x direction 
and denote the remaining d — \ directions by r . 

At the Burgers's fixed point, i.e., for large ft, x scales 
as r\_ , and p scales as r" , where a = 1/2 and z = 3/2 for 
the 2 + 1 dimensional realization of the model j|] . From 
Eq. ( [4.2[ ) one sees that the terms proportional to and 
[dxPY scale like r"^^^ and r^"~^^ respectively, which 
implies that they are irre leva nt on large length scales. 
Thus, the eigenstates of {L2) satisfy the 1-1-1 Burgers 
type equation, where the first order derivative with re- 
spect to time is multiplied by ft. For more details see 



This picture is e ssentially valid in higher spatial di- 
mensions. Eq. (4.2) is not specific to any dimension, so 
it gives the correct scaling behavior in terms of a and 
z. The fact that a + 2 = 2 in any dimension, where a 
decreases with the dimensionality of the system, implies 
that the irrelevant terms become even more irrelevant 
near the Burgers's fixed point for more than two spatial 
dimensions. Thus, one should expect that the delocal- 
ized ground state in d -I- 1 spatial dimensions has the 
scaling properties of a directed polymer in d dimensions 
while extended in the direction of the tilt. For d > 3 
this implies that for weak disorder the system may be in 
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t he w eak coupling regime, so that the nonhnear terms in 
(4.2) are both irrelevant, and the system flows into the 
Gaussian (Edwards- Wilkinson) fixed point. As a result, 
the wandering exponent of the low-lying states becomes 
1/2, like a classic lattice random walk, and the energy 
barriers between these eigenstates do not scale with the 
size of the system. However, above a critical threshold 
value of the disorder strength, the system renormalizes 
into the strong coupling regime, and a scenario similar to 
the 2 + 1 dimension case should apply. 



B. Energy statistics and finite size effects in the 
extended phase 



In this subsection we address the interpolation be- 
tween the localized and the Burgers, directed polymer 
like, limit. As the tilt h becomes strong, there is a 
linear response of the bulk magnetization to the tilt; 
the average perpendicular component of the magnetiza- 
tion bj^ = B^0o/47r is related to the external field as 
= (T(hx)h^ where cr{h±) ^ as he, and 

CT(h^) ^ 1 as hj_ — > oo. 

For hj_ > h^, the vortex line proliferates superkinks 
of the same sign II . Thus, the energy associated with 
the the existence of the flux line at point r is not U{r), 
but rather f/(r) = [/(r)r(r, J7(r)), where r(r, ?7(r)) is 
the longitudinal length, the flux line spent on the point 
r, which is itself a function of the potential energy at the 
point. Since the flux line tends to spend more time on 
strongly attractive points, U{r) is a nonlinear function of 
U{r) , and the energy distribution function P{U) is no 
longer symmetric; instead, it tends to emphasize the at- 
tractive regions while deemphasizing the repulsive ones. 
Calculation of P{U) in some limits of a simple model are 
presented in Appendix B. It is shown that the effective 
disorder strength diverges as {h~hc)^^^^ as h approaches 
the critical field from above. 

Since actual systems are of finite size, it may lead to 
completely different response to the external drift. In the 
directed polymer phase, the disorder strength and the tilt 
modulus lead to a length scale Lc, above which the poly- 
mer adapts itself to the pinning potential. Below Lc, the 
line remains stiff, since the fiuctuations of the pinning 
energy grows only sublinearly with length [Q. Thus, in 
the Burgers's phase, the low lying eigenstates may be 
fluctuating, depending on the actual length of the sam- 
ple and the strength of the disorder. As h ^ he from 
above, the effective disorder strength diverges, such that 
the ground state is always fluctuating in this limit. The 
various crossover regions for different sized samples are 
illustrated in Fig. 2. 

It is interesting to note that the critical length Lc may 
be estimated by two different ways [|l4| . First, one may 
extract the relevant part of the Eq. ( |4.2D (where we omit 
the bias b, since it does not effect the results) 




FIG. 2. "Phase diagram" for the different states of the 
flux line. In order to see the localized phase, the longi- 
tudinal dimension of the sample should be longer than 
the localization length, which diverges like {he ~ h)~^ as 
h ^ he from below. Above the transition, the potential 
barriers between two competing metastable states scales 
with the length of the sample only if L|| > Dh/A^l^. As 
h — > oo, L|j saturates to its value for point disorder in the 
xy plane. When the convection is smaller, the effective 
disorder becomes stronger; as a result, Ljj ~ (/i — hc)^^^ 
near the transition. 



h^^Dj_Vl^+-iV^cf>f+^V{r^,^). (4.3) 
ox 7 

Division by h and rescaling of the length units by some 
short length cutoff Iq, i.e., y = x/lo, one gets the dimen- 
sionless form, 



^^l<f>+^{^±4>r + ^V{y^,y^^). (4.4) 
7/1/0 



d(j) _ D 
dy\\ lo h 

The Fourier transform of 0(y±,y||), (j)(k,uj) satisfy the 
integral equation, 

( k, w) = G'o(k, Lu)V{k, uj) (4.5) 
D f dqdn 



-fhl 



- / ^ [q • (k - q)] <A(q, ^) <^(k - q, c - f7) 



where 



Go(k,c^) 



iu - 4f k2 ' 

7/1/0 



(4.6) 



and V{\i,LL!) is the Fourier transform of the growth rate. 
A loop renormalization group correction to the dimen- 

sionless diffusion constant D = — S- comes from the self 

7/1(0 

energy diagram and is of order J dgq"^^^ where the di- 
mensionless coupling constant g — -jyj^- In the (2 + 1) 
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dimensional realization of our system, Eq. (4.4) is re- 
duced to 1 + 1 dimensions, and this correction diverges 
like / q as the incoming momentum vanishes; this im- 
plies that the transverse dimension of the system, L_l, 
should satisfy Lj_ > -^^k in order to sec the effect of the 
disorder, i.e., the transition from stiff ground state to a 
wandering one. A trivial scaling argument then shows 
that the longitudinal dimension of the system is given by 
L,,|| = {Dh/^Hlf/\ 

An alternative way to see the same result is the follow- 
ing: consider a 2D isotropic homogeneous tight binding 
model, which is a discretization of Eq. (2.3). The hop- 
ping term, thus, is D/Iq, where Iq the lattice constant. 
The energy spectrum of this model is simply E{kx, ky) = 
[cos(fc2;/o) + cos(/cj,/o)], where ^ < k^^y < tt. Thus, the 
minimal separation between nondegenerate energy levels 
is AE ~ -p-. Introducing the tilt h/D pointing in the x 



direction, leads to a simple shift of k^ 



'D 



, and the 



corresponding states are tilted in the x direction, with- 
out any wandering in the direction transverse to x. As a 
result, the eigenstates are now complex, and the minimal 
imaginary distance between energy levels is P% . Now 

let us introduce disorder into the system. Using first 
order perturbation theory one finds that the leading or- 
der correction to the energy levels, is proportional to A^; 
thus, we have the same result, i.e., that non trivial effects 
of the disorder may be seen if L± > L±^c ^ {Dh/ A'^Iq). 
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APPENDIX A: TAIL STATES STATISTICS 

In this appendix we calculate the density of states 
(DOS) at the tails of the band, i.e., near the ground 
state. We then use the results to study the decay of 
these states to the ground state. The model we use is 
a discretization of Eq. (2.3), where we neglect the con- 
stant bias of the growth rate b, since it does not effect 
the statistics of the DOS. The discussion below assumes 
localization of the wavefunction of the relevant part of 
the eigenvalue spectrum, hence, the tilt term h does not 
change the spectrum; it effects the decay rate via the 
modified localization length, as discussed in Section 3b. 

We consider, thus, a model of a c? dimensional config- 
uration space tilted by square blocks of length Iq, where 



the potential energy in each of these blocks is a con- 
stant number taken from a square distribution between 
and A (the choice of these limits is only for the sake of 
convenience, since the results are independent of the po- 
tential bias 6, on the other hand, other probability mea- 
sures may give different results. In particular, for any 
unbounded distribution, such as Gaussian distributions, 
the spectrum of 7i is unbounded from above, and there is 
no delocalization transition in the sense we discuss here). 
The dimensionless Hamiltonian is : 



U(x),. 



(Al) 



It is easy to see that the DOS function p(e) in this model 
is bounded from below by zero, and p(e) — > as e — > 0. 
The tail of the density of states is determined by the 
range of parameters in which rare events, such as large 
spatial regions with low potential energy, determine the 
DOS. 

Let us estimate these fluctuations in the following way: 
the probability to find a hypersphere of radius R which 
contains only blocks of potential energy less then Vq is 

P(i?, U < Uo) - exp([i?//o]''ln[Uo/A]). (A2) 

The ground state energy of these fluctuations is given 
approximately by 



D 

ml 



(A3) 



so that the probability to get the energy between e and 
e + de using a sphere of radius R is. 



(A4) 



This expression is well defined for ^/D/e < R < oo. Op- 
timizing p(-R, e) with respect to R gives, up to logarithmic 
corrections, a maximum at the lower limit R* — i/D/e, 
so that as e ^ from below, 



(A5) 



We assume that the DOS is proportional to P{e), i.e., 
at the tail of the distribution we have g(e) ~ gQP{e), 
where go is some normalization proportional to the DOS 
at the middle of the band. 

Let us consider now the tight binding analog of the 
above model. The on site potential is now Wi, taken 
from a square distribution between zero and A. The 
Hamiltonian is 



^ = 9Im("x+e„ax + /^.c.)- 



^ Vial-G:^ (A6) 



where a^, Ox are boson creation and annihilation opera- 
tors, are unit lattice vectors, and the hopping element 
t is normalized by D. 

The eigenenergies of this Hamiltonian are bounded, 
— t < e„ < A + 1. The states at the tails correspond to a 
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rare spatial fluctuations of w^'s; in general, the chance to 
find such fluctuations (e.g., a region of radius R in which 
the on site potential is less then V) is the same as in the 
previous model. The energy spectrum of such fluctuation 
is given approximately by 



e = V + t^cos(/c^) 



(A7) 



where kr, 



thus, the tail states obey the relation 
e^^V -t + t/R^ (A8) 



so that the result (A5) is applicable here also, with the 
energy measured from the lower bound of the band. 



APPENDIX B: MAGNETIC RESPONSE AT 
LARGE TILT 

In this Appendix we study the response of the bulk 
magnetization B± to the tilt of the external field. In the 
localized phase there is no linear response, i.e., no bulk 
tilt of the vortex as a result of the external perpendicular 
magnetic field. Thus, our discussion here is relevant only 
in the extended phase. 

It is important to note that for an unbounded distribu- 
tion of defects strength, such as a Gaussian distribution, 
one flux line is always localized. For a given tilt, there is 
always a strong enough defect such that the localization 
length is small enough and the tilted state is localized. 
A delocalization of one flux line may occur only if the 
distribution of the defects is bounded. We will take the 
quenched random potential from a distribution of width 
A, so that the minimal localization length is of order 
1 / \/A. The tilt is assumed to be strong such that all the 
states are localized. 

The typical situation, though, is that the flux line 
"lives", for a while, on defects which are quite strong 
with respect to their neighbors, i.e., there are no stronger 
defects in the vicinity of the strong one. After a while, 
however, the flux line is going to jump to another de- 
fect which is weaker, in order to satisfy the tilt term in 
the free energy. This process may be described by the 
optimization of the free energy of the "local jump" 



<5F= - 



(1 + AgOi)R'i) 



- / H.-B{z)dz, 



(Bl) 



where r is the length tha t th e flux line spend on the 
defect. The first term in (Bl) indicates that there is a 
decrease of the free energy due to the fact that the flux 
line lives on the defect. This decrease is proportional to 
the potential energy difference between the current site 
and the site the flux line jumps to; it is bounded from 
above by A while it decreases to zero when the jump dis- 
tance, R, is much bigger than the density of states per 
unit area, since then the flux line may find another strong 



defect. The denominator of the first term in (Bl) reflects 
the simplest interpolation between these two limits. The 
term H • B indicates the free energy per unit length de- 
crease since the flux line average direction B tends to 
be in the direction of the external magnetic field H. In 
terms of the "local jump" , this term takes the form 



H . B = ^{r' + R'Y'^ cos{atg[^] - atg[?^]}. 



An 



(B2) 



Now let us calculate the effective tilt of the fiux line, 
i.e., the response of the internal magnetic field B to the 
external field H. For H in the direction of the defects, 
B is parallel to H; this is also true if H is perpendicular 
to the defects. Small tilt of H from the defects direc- 
tion is not enough to tilt the flux line, so that there is 
no response of the internal magnetic field (except near 
the surface, see [^). Above the critical tilt there is linear 
response of the flux line, i.e., B± = (j{H±)H±, where 
(j{H±) ^ as H± Hj^ and <t{H±) approaches 1 as 
H± — > oo. Using our local jump model, near H^^ i.e, 
where H± << and R << t; near the critical tilt 
i? ^ 0, so that we may assume that R'^ << Ag. SF then 
takes the form 



6F- 



-At- 



H4>o ,Hj^ R 2 



8tt 



Optimizing with respect to r one gets 

A. 



(-n 

T 



(B3) 



(B4) 



Since, in this limit B± ~ R/t and H± 
using the definition ei = -^1^, we have 



9 9 A 

B\ = Hl 

ei 



Hj_/H\i, and 



(B5) 



i.e., below HI = \/ A/e there is no response to the ex- 
ternal tilt while above HI B± ^ {H± — H^)^ where 
(3 — 1/2. This analysis is valid for the region in which 
both B± and H± are much smaller than one, and the 
resulting R is small compared to {gAy^'^. 

In the opposite limit, where the external magnetic field 
is perpendicular to the defects, R is much bigger than 
both T and (gA)^/"*, thus the first term is negligible and 
the fiux line is parallel to the external field. 

The above treatment was for one jump, i.e., it ne- 
glects global effects of the defects. For the actual B±, 
one should take the typical defect instead of the extreme 
one. However, we may use the same analysis in order to 
get the "time" , r, spent by the flux line on a particular 
defect, by dividing the radius of the defect by the "local 
velocity", i.e., the associated perpendicular part of the 
internal magnetic field, associated with it. The result is 
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^' /^2 



i?0 



(B6) 



where i?o is the radius of the defect. Thus, the effective 
energy associated with a defect of potential energy U is 
tilt dependent; as one approaches the critical tilt for this 
defect, H± Hj^+S, the effective energy Eeff = Ut{U) 
diverges as l/VS. (This is in agreement with the result 
for he from the localized limit, i.e., by equating the typ- 
ical inverse localization length \J~KJD to the hjD term 
in the gauge transformation |15[). 
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